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Abstract. It now becomes apparent that condensed DNA toroid which emerges 
in a poor solvent condition can be realised in the framework of the non-linear sigma 
model on a line segment. In fact, the classical solutions of the model, i.e., of the 
bending potential exhibit toroidal forms, and fit well with the delta-function term 
of the attractive interaction. The proposed theory is in good agreement with 
experimental observations that the toroid is the ground state. In this paper, 
we give a rigorous proof that the solutions are indeed the exact solutions of the 
equations of motion with the first derivatives of the attractive interaction term. 
We also show a refined mapping to experimental data, considering the finite size 
effects of the cross section and of the chain length. 



PACS numbers: 87.14.Gg, 87.10.+e, 64.70.Nd, 82.35.Lr 
1. Introduction 

Since in living cells DNA is often packed tightly, for instance, inside phage capsids, 
DNA condensation has drawn much attention [1-7] as well as its biochemical/medical 
importance in the emerging field of gene therapy [6,7]. In fact, when we put condensing 
agents as multivalent cations into DNA solution, it can cause DNA to undergo the 
condensation from a worm- like chain (whip) to toroidal states [1-7] . 

A double stranded DNA chain can be modelled by a semiflexible homopolymer 
chain [8-10]. To increase our understanding of such toroid condensation and "whip- 
toroid transition" , semiflexible homopolymers in a poor solvent condition (i.e effective 
interactions between polymer segments are attractive) have been investigated as simple 
models [11-26]. Simulations using Monte Carlo, Langevin approaches or Gaussian 
variational method, calculated phase diagram for the semiflexible chain in a poor 
solvent [11-17]. In most analytic works, existing phenomenological models balance the 
bending and surface free energies to estimate toroidal properties [18-26]. It becomes 
probable that toroid is the stable lowest energy state — the ground state. We note, 
however, that the analytic aspects of the works assumed a priori toroidal geometry 
as the stable lowest energy state with no analytic proof [23]. Moreover, compared 
to the theory of coil-globule transition of flexible chains [8,27-30], which are well 
described by field theoretical formalism [27, 28], there was no simple microscopic 
theory, which contains the salient physics to demonstrate the whip-toroid transition 
of the semiflexible polymer. 

In our previous work [31], we introduced a novel microscopic model for the 
toroidal condensation and the related conformational transitions of the semiflexible 
homopolymer chain, using the path integral method in the framework of the non-linear 
sigma model on a line segment. The proposed theory resolved the difficulties arising 
from the local inextensibility constraint of the polymer chain [9, 10] and from the non- 
local nature of the attractive interaction along the chain. The classical solutions of the 
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non-linear sigma model fit well with the attractive interaction term between polymer 
segments, and demonstrate the whip-toroid transition of the semiflexible polymer 
chain. Moreover, predicted mean toroidal radii are in quantitative agreement with 
experimental data of DNA toroids. However, it has not been proved that the classical 
solutions are the exact solutions of the equations of motion of the whole system, i.e., 
of the sigma model with the non-local attractive interaction term. 

In this paper, reviewing [31], we investigate a single semiflexible homopolymcr 
chain with attraction, using the path integral method and the nonlinear sigma model 
on the line segment. We provide a rigorous proof that the solutions of the model are 
indeed the exact solutions of the equations of motion of the system, and renew the 
derivation of the attractive potential energy. A new and refined mapping of the mean 
toroid radii to experimental data of DNA toroid is also given, considering the finite 
size effects of the cross section and of the chain length. Note that a detailed discussion 
on the stability of the solutions is given in [31], so we are not going to discuss this 
issue in what follows. 



2. Preliminaries 



In the continuum limit, the Green function (end-to-end distribution) of a semiflexible 
polymer chain with attractive interactions can be given by the path integral: 

pr(L)— R,u{L)— Uf 

G(0,R;u l ,u f ;L,W)=Af- 1 y[r{s)]e~ n ^ w \ (1) 

Jf(0)=0,u(0)=ui 

with a local inextensibility constraint |u(s)| 2 = 1 [9,10]. s is the proper time along 
the polymer chain of total contour length L. f(s) denotes the pointing vector at 
the 'time' s in our three dimensional space while il(s) = corresponds to the 

unit tangent (or bond) vector at s. Af is the normalisation constant. Following 
Freed et al. and Kleinert [9,28], the dimensionless Hamiltonian can be written 
by H [r , m, W] = J Q ds [H(s) + Vat(s)] where H(s) and Vat(s) are the local free 
Hamiltonian and the non-local attractive interaction term, respectively: 

(2) 



H(s) = l - 



dS u{s) 



Vat(s) = -W f ds'6 (\r(s) - f(s')\). (3) 
Jo 

Note that the interaction term is non-local with respect to the proper time s of the 
path integral. The persistence length / is assumed to be large enough to realise its 
stiffness and W is a positive coupling constant of the attractive interaction between 
polymer segments. The precise definition of Vat will be given at the beginning of the 
following section. Thermodynamic (3 = l/(fc_BT) is implicitly included in I and W. In 
what follows, we express r by the unit tangent vector u and therefore the Hamiltonian 

7i(u) in terms of u. Hence, the Green function G (q, R; Hi, uj; L, W^j becomes a path 

integral over u with the positive coupling constant W, regardless of r, 



G = V[u{s)] S (7 Q L dsu(s) - if) e~ H ^' w \ 



(4) 



where we used f(L) = J Q ds u(s) and the Jacobian is absorbed by Af which is neglected 
here. The delta function selects out the end-to-end vector. When W = 0, our free 
dimensionless Hamiltonian with the constraint |u(s)| 2 = 1 can be interpreted as the 
low energy limit of a linear sigma model on a line segment, or quantum equivalently a 
nonlinear sigma model on a line segment, rather than some constrained Hamiltonian 
system. 

We begin with 0(3) nonlinear sigma model on a line segment which is nothing but 
a quantum mechanics of a limited time s e [0, L] with the constraint. The constraint 
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\u\ 2 = 1 restricts the value of u onto a unit sphere S 2 . Its spherical coordinate 
decomposition in the free Hamiltonian H,(u) = ^ ds \du(s)\ 2 gives the action: 

S[6 U , <Pu] = ^J Q ds [{de u f + sin 2 6 u {dip u f] . (5) 
Minimizing the action in terms of 9 U and ip u yields the classical equations of motion: 

K = o, 



2 sin26» u 2 



[a 2 + 2(^ tl )cot0 ll a] =o. (6) 

We then find the following two types of classical solutions: 
u(s) = const, 
or 

9 = — and ip u = as + 6, (7) 

where a,b are constants. By symmetry argument, we proved that the above solutions 
represent the general solutions [31]. That is either a constant u(s) (rod solution) or a 
rotation at a constant speed along a great circle (toroidal solution). Note that they 
appear as minimal conformations in the three dimensional r space rather than in two 
dimensions of well-known Euler's elastica. In what follows, without loss of generality, 
we set a > and 6 = unless otherwise stated. 

Our aim is now to explore classical solutions and to study its corresponding 
conformations in the presence of the attractive interactions. Before we achieve this, 
we treat more carefully the attractive potential term ([3]) in the equations of motion. 
It might be obvious that the rod solution in (J7J) gives zero to the first derivatives of 
the attractive potential, because any small change of the solution contributes nothing 
to the potential. Also, the second solution seems to give zero to the first derivatives 
of the potential. Because, since its corresponding conformation sits on a single circle, 
the solution would locally give a bottom of the potential. In the following section, we 
explicitly prove such zeros of the derivatives in a rigorous way. 



3. The attractive potential part of the equations of motion 

In this section, we show explicitly the first derivatives of the attractive potential which 
enter the classical equations of motion. We then verify that they vanish in the classical 
configurations ([7]) and therefore that the configurations are in fact the exact solutions 
of the equations of motion of the whole system. The precise definition of the attractive 
potential term is given in [31] by: 



/ dsV AT {s) = -W I ds / ds' S re (R(s,s')) , 

JO Je JO 



where 



R(s,s') = \r(s)-f(s')\ = 



dt u(t) 



(8) 



(9) 



e is the ultraviolet (short range) cut-off and S re is the renormalised delta function 
which matches Kronecker's delta in its discrete correspondence. In our case, it is 

S re (R(s, s')) = R'(s, a') S (R(s, a')) , (10) 

where R' = -frR. 

ds' 

The first 9 derivative of the attractive potential can be given by: 
S ' L 

se(a ) 

= / da H step (a - s ) ds' H step {so - s)— — -S re (R(s, a')) 

Je Jo Ot){S ) 



ds I ds' S re (R(s,s')) 

s 
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ds 



max(sQ 1 e) 



7mn(so,s- 



ds' 



56(s 



-5 re (R(s,s')), 



(11) 



where H st e P (s) is the Heaviside step function. The first derivative of the renormalised 



delta function (|T0|) can be further reduced to 
Substituting (fT5|) into (TlTj) , one obtains 



S(R(s,s')) 



SR(s,s') 
S9(s ) 



(12) 



/ dsV AT (s) 
Jo 




S(R(s,s )) 
6(R(s,0)) 



max(so,e) 



5(R(s,s-e)) 



( 6R(s,3 Q y 
V 60{s o ) 

sR{ s ,oy 

S6(s ) 

( 5R(s, s — e) 



V S0(s ) 



(13) 



Similarly, the first p derivative of the potential can be given by replacing all the 9(s ) 
derivatives with those of ip(so). It should be stressed here that when the chain has 
no interactions between segments, i.e., 5(R) — 0, all the above terms trivially vanish. 
Therefore, the equations of motion are given by eq.©. For the detailed derivation of 
eqs. ([T2"| IT3|) . and another expression of (fl3|) . refer to |Appendix A| 

First, let us calculate the first 9 derivative of the attractive potential in the 
classical configurations. A naive calculation of the first derivative of R shows that 
it is, in general, of the form 1/R and becomes singular for nonzero values of 5(R). 
Therefore, it sometimes fails to estimate exactly the first derivative of the potential. 
This is indeed the case for the 9 derivative. Instead, we take rather an intuitive way 
for the exact estimate of the 9 derivative. 

Around the classical configurations, 9 is always around 7r/2, thus the small 9 
fluctuation at sq always results in an infinitesimal shift of R(s\, S2) in the 9 direction 
for Si > sq > S2- Especially when R(s\,S2) — 0, the fluctuation can be expressed by 
the sinusoidal perturbation 

6R( Sl , S2 ) = \sm(e(s )-^)\, (14) 

which clearly shows that the function R is not smooth at R = 0. Substituting 
9(s ) = 7t/2 + S9(s Q ) into fig]), we have 

5R( S1 ,S2) = |sm(<50(s o ))| = \S9(s )\. 

Therefore, 



(15) 



SR(si, s 2 ) 
59{s ) 



sgn (9(s ) - 9 cl 

assical (so)), 



(16) 



;(so) = 7r/2. By the standard definition of the sgn function, sgn(0) = 0, 



it follows that 



SR(s 1 ,s 2 ) 



S9(s ) 



(17) 



0(s o )= 



Hence, the first 9 derivative vanishes in the classical configurations: 



S9(s ) 



dsV AT (s) = 0. (18) 

classical 

An infinitesimal shift in the ip direction to R is not as simple as the one in the 
9 direction. Rather, a straightforward calculation shows that the first ip derivative of 
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the potential vanishes as well in the classical configurations. By substituting |(7J) into 
([5]), the function R can be explicitly written down for the classical solutions 0: 

2 a, 
- sin -(s-s ) . 



R{s,s') 



classical 



Similarly, its tp derivative is given by: 



SR(s,s') 



8ip(s ) 



cos a(s — so) — cos a(s' — s ) 



(19) 



(20) 



classical 

When R approaches zero, that is, s = 
the l.h.s. of eq.(|2U)) becomes 

cos a(s' + £ — Sq) — cos a(s' — Sq) 



2|sinf(s-s')| 

- 22™ + £ w ith \e\ < 1 and m e Z+, then 



l.h.s. 



2 | sin f | 



f/5 



cos a(s —sq) + sgn ysm — J cos 



(75 



■s ). (21) 



Therefore, the first tp derivative turns out to be well-defined. For the specific values 
of s' , one finds 



8R(s, s ) 
Sf(so) 

SR(s,0) 
8tp(s ) 



ae 

Y 

as 



sm ■ 



cos(aso) — sgn ^sin 



ae 

cos — sm(aso). 



(22) 



Now we are ready to calculate the derivative around a singularity of 5(R(s, sq)): 



de S 



as 



sm ■ 



(IS 



sm ■ 



dX S(X) 



l\x\ 



1 - 



i, 2 X' 2 



= 0, (23) 

where X = 2gi n (ae/2), 6 and 8' pick up the neighbourhood of the zero only (i.e., 
< 8 <C 1), and the identity #(|X|) = 8(X) is used. The same holds for the first term 
of R(s, 0) in (|22p. As for the second term, one should be aware from the existence of the 
sgn function that it behaves odd around the zeros of the delta function. Accordingly, 



this term does not contribute to the final result. Finally, i?(s, s — e) , , = e > 

^ ' v 7 ! \ classical 

and 8(e) = so that the third term in (|13jl vanishes. Hence, the first tp derivative of 
the attractive potential also vanishes: 

i-L 

dsV AT {s) =0. (24) 

classical 



8tp(s ) Jo 



4. Whip, toroid, and the potential energy 



Now let us consider the attractive potential ^ itself. It is difficult to interpret it in 
the context of quantum theory due to its non-local nature along the polymer chain. 
However, we can derive its explicit form in the classical configurations l[7]). Substituting 
l[T9"]) into Ijgj), one finds 

r-L 

ds V AT {s) = - W I ds I ds' 8 re 



W 
W 

w 



( sin I (s - s') 



ds 



ds 



ds' 8 re 



o 

as 
2^ 



(s — s ) mod — 
a J 



= - w 



N(L)-1 

-- > : k+- 

a ^ 

k=l 



2n 
a 
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= —N(L) (ML) + 1)-WL N(L), (25) 
a 

where 5 re (const x x) — 8 re (x) is used, and N(s) = [as/2?:] by Gauss' syniboQ. 
ae/2iT 1 is assumed without loss of generality. Note that N(L) represents the 
winding number of the classical solution ([7]) along a great circle of S* 2 . Therefore, 
the form of the potential can be thought of as 'topological' terms from the winding 
number. 

The non-zero winding number of the classical solution in the u space means that 
the polymer chain winds in the r space as well. That is, when a > configurations 
around the second classical solution start forming a toroidal shape since 
/ I {sin(as + b) — sin(6)} \ 
f(s) = -i{cos(as + 6) -cos(6)} , (26) 
V const. J 

and stabilise itself by attracting neighbouring segments. We call such classical 
solutions the "toroid states." Whenever a increases and passes through the point 
^jp- for n S Z_|_, another toroid state appears with the increased winding number n. 
Note that the radius of the toroid state is given by -. When the chain of contour 
length L winds N(L) times we have the N(L) circles of each length and the rest 
\L — ^-N(Vj). The first and second terms in eq. (|2"5"|) result from the former and the 
latter, respectively. When < a < the chain cannot wind like the toroid states. 
Both ends of the chain are not connected to each other, thus can move freely as well as 
any other parts of the chain fluctuate. As long as the total energy of the chain does not 
exceed the bending energy of ^j^- at a = they can whip with zero winding number. 
We call such low-energy states the "whip states." Although the definition includes 
fluctuations around the classical solutions, unless otherwise stated, we primarily refer 
to the classical solutions of such states, which are rather bowstrings than whips. 



Whip states 



Toroid states 





N = 2 



N=0 



Figure 1. The whip (N = 0) and toroid states (TV > 1). 



5. Energy levels of the analytic solutions 

We now explore the exact energy levels of the whip and toroid states, and briefly 
mention the conformational transitions between these states. The dimensionless 
Hamiltonian of the classical solutions ([7]) is a function of l,L,W and a: 

H c i(a,l,L,W) = ^-a 2 + —N(L) (N(L) + l)-WLN(L), (27) 
2 a 

with the definition of N ^ = for a = 0. It is convenient to rewrite this in terms of 

a 

the conformation parameter c = -K (^) and x = lp, and thus with [x] = N(L): 

2tt 2 / 

H d (a,l,L,W) = WLH{c 7 x) = —— (4cH(c,x)), (28) 

1j 

| Gauss' symbol [x] gives the greatest integer that is not exceeding x. 
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where 

«M^ + ^-H. (29) 

Therefore, the function H(c, x) can solely express the shape of the energy levels. 
Without loss of generality, we assume x > 0. 

Let us consider a case with L, W, and / fixed. By definition, 7i(c, x) is continuous 
in the entire region of x > and is smooth in each segment: x € [n, n+1] for rt 6 Z>o- 
However, it is not smooth at each joint of the segments: x € Z+. We plot, in FigO the 
energy levels as a function of x for different values of c, showing qualitative agreement 
with Conwell et al. for DNA condensation [6]. Given [a;] = N(L) is fixed, the 



Energy levels 























































c= 1/2 










(=27/16 
c = 9/4 
c = 4 





















1 2 3 4 5 



Figure 2. The dependence of the energy AcH(c,x) on x = aL/2-n and c. The 
unit in energy is set to be the maximum bending energy of the whip state. 

Hamiltonian ([HD takes a minimum at x = x c (N(L)) = (cN(L)(N(L) + 1)) 1/3 . When 
N(L) < x c (N(L)) < N(L) + 1, TL(c,x) behaves quadratic in x and x = x c (N(L)) 
indeed gives its minimum in the iV(L)th segment. In other words, the condition of a 
(meta-) stable state in the iV(L)th segment turns out to be 

N L (c) < N(L) < Nu(c) for c > 4, 

1 < N(L) < N v {c) for < c < 4, (30) 

where N L (c) = f (l - § + , and Nu(c) = f (l + ^/l+§) . 

When Nu(c) < 1 (i.e., c < i), the above second condition vanishes and thus the 
whip states only survive at low energy. Therefore, at the critical value of c = ^, the 
whip phase to whip-toroid co-existence phase transition may occur. When c > i, there 
always exists at least one (meta-) stable state. The number of minima are roughly 
given by the width of the region for N, i.e., Njj(c) — Nl(c). For example, when c > 9/4, 
either Nu{c) — Nl{c) > 3 or Njj(c) > 3 with the lower bound being 1, hence there 
are at least three minima with different positive winding numbers. This occurrence 
of the multiple local minima and stable states is because the energy is given by the 
balance between the bending energy and the attractive potential energy. The former 
is a monotonically increasing function of x while the latter is monotonically decreasing 
but not smooth function of x. Note that the number of minima may be reduced in 
some cases, for example, with some finite size effect on. The existence of multiple 
minima indicates the first order transitions between these stable states. In fact, when 
the energy of the N = 1 stable toroid state becomes zero, the whip dominant to toroid 
dominant phase transition occurs. Such a value of c is 27/16. When c > 27/16, the 
toroid states will dominate the action. Further discussions on the transitions can be 
seen by changing the value of c and are the subject of a forthcoming publication. 
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Similarly, one can deduce the inequality for the critical values of c: 

cf ) <c<cW (31) 

where c^= and c[j = ^-jp—- In this region, the iVth segment has a minimal 
and (meta-) stable point. From this, one can approximate the relation between the 
conformation parameter c and the dominant toroid winding number N c . Furthermore, 
it leads to the asymptotic relation for either large c or large N c : 

ci^+c^ 1 /l 



2 = iVc+ 2 +0 UJ" 7Vc - (32) 

The first (~) is the approximation and the second one denotes the asymptotic relation. 

6. Comparison to experiments and concluding remarks 

For large c, the ground state — the dominant toroid state of the winding number N c 
can be estimated by the inequality relation (f3Tj) of c. It reads N c ~ c. Using this, we 



can estimate the radius of our ideal toroid: r c — „ \ T = 4§r- Our ideal toroids have 
zero thickness, but the finite size effect of their cross sections can be approximated 
by the hexagonally arranged DNA chains with a van der Waals type interaction, i.e., 
with the effective nearest neighbour interactions. In the case of N(L) > 4, it leads to 
the modified Hamiltonian [31]: 

H cl (a, I, L, W) = ^-a 2 - ^V(N(L)) 
2 a 

a \ 27r J 

where V(N) = 3N ~ 2y/Zy/N - 1/4 and Gap(N) = V(iV + 1) - V(iV). Following the 

same procedure for < c < , we find iV c ~ (2-\/3c) 5 for large c. By defining 
the mean toroid radius (i.e., the average of inner and outer radii) by r c = we 
now estimate the mean radius of the toroid in a physical system in detail. A coupling 
constant can be given by W — where k is the number of the electric dipoles 

in a monomer segment, which create van der Waals interactions of the magnitude 
e w . l m denotes the length of the monomer along the chain contour, taken to be 
l rn ~ 5 bp = 1.66 nm in the end. Substituting N c and the above, we obtain the 
asymptotic relation up to numerical factor: 

r c ^(6n)-Hi^y =(6 7 rr*i*(U) f (^;) (34) 

This L~> dependence agrees with the proposed exponent in the asymptotic limit 
[19,21,24,26,32]. 

We estimate the mean toroidal radius of T4 DNA in low ionic conditions reported 
in [5]. Substituting L = 56AA/im, I ~ 50^ 60 nm, and l m into Q34p , the mean radius 
of the toroid is found to be r c = 29.03.B-i - 31.23^-i [nm], where B = This 
result is in good agreement with the experiment r c ~ 28.5 nm for B» 1.15. The same 
argument for Sperm DNA packaged by protamines (L = 20.4/xm) [4] gives the analytic 
value r c = 23.69B~s - 25.485' " s [nm] , which also agrees with an experimental result 
r c ~26.25 nm for B k, 0.851. Note that the latter has a larger diameter of the segment 
and is expected to have the weaker interaction with smaller _B« 0.851. 

The above mappings are given on the basis of the asymptotic relation Q34p . but 
we can find the range of B numerically without using the asymptotics. The relation 
r c = L/(2irN c ) means iV c is given by the inverse of r c up to constant for a fixed 

contour length L. So, the approximation, c ~ (c\f c ^ + c^"^ /2, leads to l/B as a 
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function of r r 



I 

B 



V3L 2 

7T 2 / 

/( I'm 



rc /2 {(l + r c ) v /l-r c /4- + 3f c /4} 
1 + (1 + f c ) 3 



(35) 



where r c = 2nr c /L. Therefore, we can plot the r c - / relations for the fixed values of 
B and L (Fig©. It reads B = 1.05 - 1.26 w 1.16 {B asympt « 1.15) for the T4 DNA 
and B = 0.784 - 0.940 w 0.862 (B asymvt w 0.851) for the Sperm DNA, respectively. 

r - / relation & allowed regions of 5 




Figure 3. The relations between r c and Z for T4 DNA with B = 1.053 and 
B = 1.263, and for Sperm DNA with B = 0.7837 and B = 0.9404. The two 
regions are given when B is changed continuously from one to the other. 

The exponents v predicted in the literature for r c ~ L v are v = i in most 
cases [19,21,24,26,32]. However, they are inconsistent with the experimentally well 
known observation that the radius is independent of the chain length [2-5]. This 
might suggest that the real interaction is not van der Waals like, or at least is not 
a single van der Waals type interaction. It should be noted here that combinations 
of our ideal toroid and its finite size effect can give a range of v = — 1 ~ ~ in some 
region. An interesting remark is that when we apply Coulomb like interactions to our 
approximation, we observe that the radius nearly remains constant as L changes. The 
precise analysis is to be studied and given elsewhere. 

So far we have dealt with the classical solutions, which are derived from the first 
derivative of the action. Thus, they may correspond to the global/local minima of 
the action in the configuration space. In fact, when c > 4, the toroid states become 
stable under the fluctuations, since deviations from such toroid states will cost a large 
penalty in energy. Its detailed discussion is given in [31]. Note however that, when 
c is small, a number of the whip states become equally or more probable than any 
toroid state. Another interesting question will be how the derivatives of the potential 
behave when the configuration be deviated a little, but not infinitcsimally, from the 
classical configurations. To answer this question is beyond the scope of this paper. 

Acknowledgments 

The authors are grateful to K. Binder and his group for their stimulating 
discussions and hospitality in Mainz. One of the authors (Y.I.) is grateful to K. 
Nagayama for his discussions and encouragement. The other (N.K.) is grateful to S. 
Stepanow and S. Trimper for their stimulating discussions and hospitality in Halle, 
and acknowledges the Deutsche Forschungsgemeinschaft for financial support. 



DNA toroid condensation as analytic solutions 10 

Appendix A. Relevant calculations for the first derivatives of the potential 

The attractive potential energy is symmetric under the exchange of s <-> s', since 
R(s, s') = R(s' , s). Therefore, one can rewrite or regard the double integral in (|8|) as 



ds 



s-e 1 

ds' = - 
2 



ds 



ds' 



L-e 



ds 



ds' 



(A.l) 



s+e 



On the other hand, the first derivative of the renormalised delta function can be 
reduced: 

o6{s ) 

' ] ' r ' " S>) -6 (R(s, s')) + R'(s, s')jJ-^S (R(s, s')) 



se(s ) 



se( Sa ) 



d 

ds 1 



f SR(s,s') \ s(nt dR(s,s') (1) n , SR(s,s') 

A ^77^ )HR(s,s))+ dg , SW(R(s,s)) saf ^ 



60(s ) 



66(s ) 



, d ( 5R{s,s') 



8 (R(s,s)) 



SR{s,s') 
S0(s Q ) 



_d_ 

ds 7 



(A.2) 



V S9(s ) 

where 5^'(x) is dS(x)/dx, and the s' and 9 derivatives are assumed to be independent. 
Substituting this into eq. fTT]) , one obtains 



S0(s Q ) 



f dsV AT (s) 
Jo 

ds J 5(R(s, s )) 



( SR(s,s ) 



min(sQ-\-e,L) 

min{sQ-\-€,L) 



max(so 1 e) 
L 



ds l5(R{s,s-e)) 



SR(s, s — e) 



S9(s ) 



I ds ( 

J max(so,e) K. 



ds \S(R(s,0)) 



f SR(s,0) 
\ 66(s ) 



(A.3) 



Here, if we revive the relation (jA.lj) . we find the symmetric expression for the first 
derivative of the potential: 



56(s ) J 



ds Vat(s) 

max(sQ — e,0) ( 

ds \s(R(s,s )) 







ds { 5(R(s,0)) 

!nax(so,e) 



SR(s,s Q ) 

se(s ) 

fSR{s,0) 



ds l5(R(s ,s)) 

max(sQ-\-e,L) I 



V S0(a o ) 

Similarly, the first tp derivative of the potential is given by 
5 



min(L — e,so) 

ds lS(R(L,s)) 



6R(s ,s) 



S9(s ) 

(6R{L,s) 
V S6(s ) 



(A.4) 



S(p(s ) 



f dsV AT {s) 
Jo 



max(so — €,0) ( 

ds \s(R(s,s )) 



o 



ds { S(R(s,0)) 

;nax{sQ 1 e) 



5R(s,s ) 



5ip(s ) 

(6R{s,0) 
\ S(p(s ) 



ds i6(R(s ,s)) 

max(sQ-\-e,L) I 



SR(s ,s) 



Tnin(L — e,so) 

ds id(R(L,s)) 



5ip(s ) 

f SR(L,s] 
V M s o) 



(A.5) 



References 



DNA toroid condensation as analytic solutions 



11 



[1] L. C. Gosule and J. A. Schellman, Nature 259, 333 (1976). 

[2] P. G. Arscott and V. A. Bloomfield, Biopolymers 30, 619 (1990). 

[3] V. A. Bloomfield, Biopolymers 31, 1471 (1991). 

[4] V. A. Bloomfield, Curr. Opinion Struct. Biol. 6, 334 (1996). 

[5] Y. Yoshikawa, K. Yoshikawa, and T. Kanbe, Langmuir 15, 4085 (1999). 

[6] C. C. Conwcll, I. D. Vilfan, and N. V. Hud, Proc. Natl. Acad. Sci. 100, 9296 (2003). 

[7] N. V. Hud and K. H. Downing, Proc. Natl. Acad. Sci. 98, 14925 (2001). 

[8] M. Doi and S. F. Edwards, The theory of polymer dynamics (Clarendon Press, Oxford, 1986). 

[9] H. Kleincrt, Path Integrals in Quantum Mechanics, Statistics, and Polymer Physics, and 
Financial Markets (World Scientific Publishing Company, 2004). 

[10] K. F. Freed, Adv. Chem. Phys. 22, 1 (1972). 

[11] A. Montesi, M. Pasquali, and F. C. MacKintosh, Phys. Rev. E 69, 021916 (2004). 

[12] I. R. Cooke and D. R. M. Williams, Physica A 339, 45 (2004). 

[13] H. Noguchi, S. Saito, S. Kidoaki, and K. Yoshikawa, Chem. Phys. Lett. 261, 527 (1996). 

[14] H. Noguchi and K. Yoshikawa, J. Chem. Phys. 109, 5070 (1998). 

[15] Y. A. Kuznctsov and E. G. Timoshenko, J. Chem. Phys. Ill, 3744 (1999). 

[16] Y. A. Kuznctsov, E. G. Timoshenko, and K. A. Dawson, J. Chem. Phys. 105, 7116 (1996). 

[17] V. A. Ivanov, W. Paul, and K. Binder, J. Chem. Phys. 109, 5659 (1998). 

[18] A. Y. Grosberg and A. R. Khokhlov, Adv. Polym. Sci. 41, 53 (1981). 

[19] M. R. Stukan, V. A. Ivanov, A. Y. Grosberg, W. Paul, and K. Binder, J. Chem. Phys. 118, 
3392 (2003). 

[20] N. V. Hud, K. H. Downing, and R. Balhorn, Proc. Natl. Acad. Sci. 92, 3581 (1995). 

[21] J. Ubbink and T. Odijk, Biophys. J. 68, 54 (1995). 

[22] J. Ubbink and T. Odijk, Europhys. Lett. 33, 353 (1996). 

[23] B. Schnurr, F. C. MacKintosh, and D. R. M. Williams, Europhys. Lett. 51, 279 (2000). 

[24] B. Schnurr, F. Gittes, and F. C. MacKintosh, Phys. Rev. E 65, 061904 (2002). 

[25] G. G. Pereira and D. R. M. Williams, Europhys. Lett. 50, 559 (2000). 

[26] I. C. B. Miller, M. Keentok, G. G. Pereira, and D. R. M. Williams, Phys. Rev. E 71, 031802 
(2005). 

[27] I. M. Lifshitz, A. Y. Grosberg, and A. R. Khokhlov, Rev. Mod. Phys. 50, 683 (1978). 

[28] A. L. Kholodcnko and K. F. Freed, J. Phys. A: Math. Gen. 17, 2703 (1984). 

[29] P. G. de Gennes, Scaling Concepts in Polymer Physics (Cornell University Press, New York, 
1979). 

[30] A. Y. Grosberg and A. R. Khokhlov, Statistical physics of macromolecules (American Institute 

of Physics, New York, 1994). 

[31] Y. Ishimoto and N. Kikuchi, J. Chem. Phys. 125, 074905 (2006). 

[32] A. G. Cherstvy, J. Phys.: Condens. Matter 17, 1363 (2005). 



